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Abstract 

We study the rate and (j> momentum distribution in semi-inclusive decays B — ► 4>X S 
induced by the quark level processes b — > <ps and b — > 4>sg, in which the gluon is radiated 
^ ■ from the internal charm quark loop or emitted from the virtual gluon of the strong penguin 

h ; 

(inner bremsstrahlung). We find B(b — > <fis) = 6.7 x 10 5 and B(b — > </>sg) = 3.8 x 10 5 . 
The momentum spectrum of <^> produced by 6 — ► 0s<? is very broad. With the cut |k^| > 2.0 
GeV, Bib — > 0s) = 6.1 x 10 -5 (where the Fermi motion of the 6-quark in the i?-meson 
is described by a Gaussian), and Bib — ► 0sg) = 1.0 x 10 -5 . Due to the special nature of 
(j), many difficulties which hindered a reliable theoretical prediction for B — > t/X s decay 
are absent in the process -B — ► </>X s . Therefore, theoretical predictions for — ► </>X s are 
relatively clean. Moreover, the clear experimental signature of the <fi is of great help. Data 
for B — > 0X S , both the branching ratio and the momentum distribution, would teach us 
about the strength of strong penguins which might be of great importance in the search 
for CP violation and for new physics at B factories. 
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1 Introduction 



Understanding of pure penguin decays of B mesons, is of utmost importance. Penguins have 
been serving as powerful probes for the Standard Model (SM) and for beyond the SM scenarios. 
Measurements of the electromagnetic penguin process b — > S7 by CLEO pj, which agrees with 
the SM, has provided very stringent constraints on new physics. The strong penguin processes 
B -> rfX e , B -> r/Kl*\ B -> <j>K& have been observed by CLEO @, BABAR |, |] and 
BELLE 0, 0. As is well known, it is quite difficult to provide theoretical estimates of exclusive 
nonleptonic decays, although there has been some progress on this topic. In particular we refer 
to the QCD factorization and pQCD approaches || developed recently. To explain the 
large yields of rj', we encounter unknown parameters like the content of 7/, mixing-angles, the 
gg — 7/ coupling and so on, which have hindered reliable theoretical predictions. There are also 
suggestions that new physics could enhance the magnitude of the strong penguin. At present 
we cannot conclude how large the window is for new physics or whether it is required to explain 
the data. The semi- inclusive B — > <pX s decay is theoretically cleaner than B — > rj'X s , since <fi 
is almost a pure ss state with mass larger than 2M^ and does not couple to two gluons. In 
addition, <fi coupling to three gluons is highly suppressed by the OZI rule. 

There is a number of studies of B — > K^(f>, using different approaches ||, 10, 11, 12 



Unfortunately the results do not converge. QCD factorization predictions |L0[ |TT| are smaller 
than those of pQCD |Tj|. Compared with exclusive decays, the theoretical predictions for 
inclusive and semi-inclusive decay rates of B mesons rest on more solid grounds. The semi- 



inclusive decay B — > cj)X s was studied a few years ago JT3[- Recently it has caught renewed 
interest [IT], [TJ[ by generalizing the QCD factorization formalism |7J] to semi-inclusive processes. 
In the present paper, we will study both the branching ratio and the <fi momentum distribution, 
hereafter denoted by [k^. We take into account the Fermi motion of the fe-quark for b — > s<fr 
and the b — > sg*g strong penguin effects arising from the inner bremsstrahlung processes in 
which a gluon is emitted from the charm loop or from the virtual gluon in b — > sg*. 

Possible large b — > sgg contribution to inclusive B meson charmless decays was discussed 
in Refs. |T5|, [16|, [T7|, [L8], [TjJ. Furthermore, while studying penguin effects, Gerard and Hou [|1^, 
|l(J found that the higher order processes b — > sgg and b — > sqq dominate over b — > sg. 



Subsequently, it was found by Simma and Wyler [17], and independently by Liu and Yao [18 
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that b — > sgg is considerably suppressed as compared to b — > sgg. Both groups found that the 
large form factor is absent when both gluons are on-shell, while if one of the two gluon 

goes off-shell, Fi(x) survives. 

In the present paper we work within the framework of an effective low energy theory with five 
active quarks which is obtained by integrating out heavy top and heavy gauge bosons ||20|| . We 
also use the QCD factorization framework to deal with the hadronic dynamics of (j) formation. 
We find that the b — > <psg contribution is very significant, B(b — > <ftsg) = 3.8 x 10 -5 . After 
imposing a momentum cut |k^| > 2.0 GeV, B(b — > <j)sg) is reduced to 1.0 x 10 -5 , which is about 
16% of the fast <fi production due to b — > <ps. To understand this "unusual" large contribution, 
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we note that the amplitude for b — > <psg is characterized by a factor Cij^a — 0.029 which 
is numerically comparable to Cq = —0.041, known to be the largest coefficient of the strong 
penguin four Fermion operators. It is interesting to note that similar large next-to-leading 
(NLO) corrections have also been found for b — > sg |I9] and B — > K*^ |[21|| . 

The remainder of the paper is organized as follows: In Sec. 2 we study the momentum 
spectrum of <fi resulting from b — > <ps by taking into account the effect of Fermi motion. Sec. 3 
contains a calculation of the b — > <psg contributions. We present and discuss our results in 
Sec. 4. Some useful formulas are given in three Appendices. 



2 B — ► (f>X s induced by b 



In the low energy effective theory of the SM, the relevant effective Hamiltonian is written as 

follows gg 

4G F 



H 



eff 



V2 



2 / 10 

v cb v c * s y: CtOi - V tb V* C g O g + £ cpi 

i=l \ j=3 



(1) 



The operators O are listed in App.A. The Wilson coefficients evaluated at scale ji = are [2^ 



Ci = 1.082, C 2 = -0.185, C 3 = 0.014, 

C 4 = -0.035, C 5 = 0.009, C 6 = -0.041, 

C 7 = -0.002/137, C 8 = 0.054/137, C 9 = -1.292/137, 



(2) 



Cio = 0.262/137, C g = -0.143. 
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Figure 1: Leading diagram for (B — ► (j)X s )2 defined in Eq|| Oi are strong penguin operators in which 
top penguin effects are embedded. 



In the naive factorization approach, b — > cf)s decay is a color- suppressed process depicted in 
Fig.l. The amplitude can be easily written as 



y/2 



Vi 



(3) 



where 



A p = -V tb V* s (a 3 + a 4 + a 5 - ^ (a 7 + a 9 + ai )) , (4) 
and a 2n -i = C 2n ~i + C 2n /N c , a 2n = C 2n + C 2n -i/N c , (n = 2, 3, 4, 5). The branching ratio is 



G\fW 



<t>""b 1 a ,2 



16tt 



m 



2 \ 2 



1 + 2- 



(5) 



Let us denote the two body contribution to B — > <f>X s by 



B[{B^<PX S ) 2 ] = B{b 



(6) 



Using = 233 MeV, m b = 4.8 GeV, and the Wolfenstein parameterization V t b = l,V ts = —AX 2 



with A = 0.817 and A = 0.2237 from Ref. [22], we have 



B[{B -> 0X S ) 2 ] =4.9 x 10" 



f naive factorization). 



(7) 



It is believed that naive factorization works very well for color-allowed processes since 
Bjorken's color transparency argument [23] applies, while for color-suppressed processes non- 



factorizable effects could play an important role. As demonstrated by Cheng and Soni |14[], 
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(a) 




(d) 




(b) 




(e) 





(f) 



Figure 2: Subleading diagrams for (B — > <j)X s )2 defined in Eq|| 



the nonfactorization effects depicted in Fig. 2 could be calculated in the BBNS QCD factor 
ization framework. The calculation of the effects of Fig. 2 was carried out in Refs. |TJ], |23 
Incorporating the nonfactorizable contributions, A p in Eq.[| is modified to 



a' p = v ub v: 



+ v cb v: 



a u + a u + a u__ {a u + a u + a U Q)+a U Qa 



1 

2 



"10 J 



10a 



(8) 



where the coefficients a?'s can be found in Ref. [24|. We subtract the contribution of gluon- 
spectator interactions (/// in [23), choose 7 = 54.8° and present the numerical results in Table 
1. From Table 1, we can see that our results agree with those by Cheng and Soni [I3|. We 
have 

B [(B -> cj)X s ) 2 } = 6.7 x 1(T 5 (QCD factorization). (9) 

This large branching ratio enhances the feasibility of measuring the strength of the strong 
penguin to test the SM by studying this decay mode at BABAR and BELLE. 

To study the momentum spectrum, we employ the ACCMM model |25j (for an earlier version 
see ||26|| ). In this model the bound state corrections to free 6-quark decays are incorporated by 
attributing to the 6-quark Fermi motion within the meson. The spectator quark is handled as 



Table 1: Numerical values for af (in units of 10 4 ) in QCD factorization and in naive factorization. 
Our results are given for 7 = 54.8°. Cheng and Soni's results \Lm are displayed for comparison. 



Our results 



Ref. M 



Naive factorization 



c,u 

a 3 ' 


76.0+27.8i 


(«3 


) 74+26i 


23 


a\ 


-375-71. 6i 


(a 4 ) 


-353-58i 


-303 


«4 


-318-15H 






-303 


c,u 

a 5 


-68.3-31.4i 




) -67-30i 


-46 


c.u 

a 7 


1.37+0.3i 


K)- 


■0.89-1. 13i 


1.16 


c.u 

CLg 


-90.8-1.4i 


(a 9 ) 


-92.9-2.8i 


-87.9 


C.U 

«10 


1.95+7.23i 




0.6+6.4i 


12.2 


C,M 

a 10a 


-0.52-1.0i 









an on-shell particle with definite mass m sp and momentum |p| = p. Consequently, the 6-quark 
is considered to be off-shell with a momentum dependent virtual mass W{p) 



W\p) =M A B + m sp - 2M B Jm% + p 



(10) 



in which energy-momentum conservation is imposed. The momentum of the fe-quark is modeled 
by a Gaussian distribution function with a free parameter pp 



4 1 

—=-3- exp 
kPf 



P_ 

Pi, 



(11) 



It is interesting to note that the parameter pp and the 6-quark average mass (rrib) were obtained 



from a fit to the photon momentum distribution in B — > X s ^ by CLEO |27[], to read pf = 410 
MeV and (m b ) = 4.690 GeV. Using these values and Eqs.[T(], [TT|, we get m sp = 298 MeV for 

Now, the momentum spectrum of resulting from the decay b — > cj)s of a 6-quark of mass 
W(p), is given by 



dTQk^p) W(p)T (W(p)) r 



d\k 4 



|k 5 J) 



(12) 



where _ provide the limits of the momentum range which results from the Lorentz boost, 
i. e. 



ki(p) 



W(b) 
6 



(E b k ±pE ) 



(13) 




0.5 1 1.5 2 2.5 3 



lig(GeV) 

Figure 3: Momentum spectrum of <j> in B — > (j)X s decay. The solid curve is for b — > 4>s including 
Fermi motion. The dotted curve is for b — > (ftsg. 



with 



Finally, 

To get the momentum spectrum of <fi from the semi-inclusive decay of the B meson, we have 
to fold in the 6-quark momentum probability distribution as given by 

dB(B -> <PX S ) «w 2 dr(|k^|,p) 

- r B / dp<f) F {p)p — -f- — , (16) 



d\k^\ Jo d\k 

where p ma x is determined by W(p) 2 > (m^ + m s ) 2 , leading to 



Pma* = (Ml + m% - (m + m s )f - ±Mlm%. (17) 



Our numerical results for the momentum spectrum are presented in Fig. 3. We can see 
the spectrum peaking near |k^,| = 2.4 GeV which corresponds to Mx s = 1.15 GeV. To suppress 
indirect <fi production, one can impose a momentum cut |k^,| > 2.0 GeV. With this cut, the 
branching ratio is 

B [(B -> 0X S ) 2 ] = 6.1 x 10" 5 (QCD factorization + Fermi motion+|k | > 2.0 GeV). (18) 
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3 Contribution of the three body decay b —> (j)sg 



The b — > (psg contribution to the momentum spectrum of (b has been studied by Deshpande 
et al. [0 for the case of gluon radiating from external quarks. They find that the effect is rather 
small, B(b —>■ 4>sg) / B(b — > <ps) « 3%. From now on we will neglect gluon emission from external 
lines, and study instead the effect of the gluon radiating from an internal quark or splitting 
off an internal off-shell gluon as shown by Fig.4. Both these processes are sometimes referred 
to as "inner bremsstrahlung" . The calculation of Fig.4. (a)-(b) is straightforward. However the 




s b 




(b) 




s b 




Figure 4: Diagrams for B [(B — ► 4>X S ) 3 ] defined in Eq.23 



calculation of Fig.4. (c) and (d) is tedious. A calculation of the charm loop for b — > sg*g has 
been carried out recently by Greub nad Liniger [113, in a detailed study of 0(a s ) corrections 
to r(6 — > sg). It gives 



Jas = IJCff, k){T A , T B } + T-Jq, k) [T , T ]. 



(19) 



'a/3 — ± otPV±i ru )\ ± ^ 5^ ± a(5\ 

q, A and a denote the momentum, color and Lorentz index for the on-shell gluon, respectively, 
while k, B and (3 stand for the same attributes for the other, off-shell, gluon. The quantities 
T i(<?, fc) and T- p (q, k) are |9| 

kr 
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32vr 2 



E(a, (3, k)Ai 5 + E(a, (3, q)Ai 6 - E((3, k, g)^Az 23 



-E(a, k, q)-^-Ai 25 - E(a, k, q)-^-Ai 2e 
q ■ k q ■ k 



(20) 



g 2 

T aP {q,k) = — 



ftg Q/3 Ai 2 + ^g a pAi 3 + Az 8 + lakpAiu + 7 Q g /3 Az 12 

L. (21) 



+ ^Ai 15 + fi^Ai 17 + ^Ai 19 + ^Az 21 



The dimensionally regularized expressions for the Ai functions are given by Greub and Lin- 
iger [|19|]. To reduce the difficulty in numerical calculations of multiple integrals, we special- 



ize their functions to d = 4 and integrate out the two Feynman parameters. For d = 4, 
E(a,/3, , ~f) = — ^60/37^7^75, in the Bjorken-Drell conventions. The analytical expressions for the 
functions Ai are collected in App.B, where we have used the MS scheme. 
The dimensionless variables s, t and u are defined as 

„ _ (P<t> + Vsf t = (gj + q? u = (q + Ps) 2 ^ 



m% ml ml 

We write the amplitude of the three body contribution to B — > <pX s as 

B[(fl-^X.) 3 ] =B(6->^). (23) 

Now, neglecting gluon bremsstrahlung from external quarks, the inner bremsstrahlung diagrams 
(see Fig. 4), lead to the matrix element 

M 3 = M 0g + M Qc + M A+ + M A _ , (24) 

and its square reads 

I -M 3 I 2 = \Mof + \M f + \M A+ \ 2 + \M A S + mMoMoJ + ^(M A+ Mi_) 

+2&(MA + Ml g ) + 2&(Ma + MI c ) + 2^{M A _M ] 0g ) + 2^{M A _M ] 0c ). (25) 

The explicit forms for the terms in the amplitude squared are given in App.C. 
The decay distribution of the <f> in B [(B — > <pX s ) 3 ] is 

dudt 256tt 3 6 2 1 1 6 1 1 ' V ; 

where the factor | is due to the average of spin and color of the 6-quark. In the phase space 
integration, we impose the cut (p s + q) 2 > m\. 
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4 Results and discussion 



Our results for the contribution of b — > (psg to the spectrum, are displayed in Fig. 3 by the 
dotted curve, which is much broader than that due to b — > <ps. One can improve the predictions 
for the spectrum of the <p mesons by adding Fermi motion to the three body processes, including 
radiation from external quarks. The contribution of b — > <j)sg is quite substantial, 

B[(B ^ ( f ) X s ) 3 } = 3.8 xlO" 5 , (27) 

where the left hand side of the last equation was defined in Eq.|2"3"|. Now, after applying a 
momentum cut |k^| > 2.0 GeV, the branching ratio is reduced to 

B [(B -> 0X,) 3 ] = 1.0 x 10" 5 > 2.0 GeV). (28) 

which is about 16% of Bib — > <ps). Interestingly, such a large NLO correction was also found 
in b — > sg [H| and B — > K*^f [^TJ. In Ref. ||19|| , Greub and Liniger have found that the next- 



to-leading lorgarithmic result B NLL (b — > sg) = (5.0 ± 1.0) x 10~ 5 is more than a factor of two 
larger than the leading logarithmic result B LL (b — > sg) = (2.2 ± 0.8) x 10~ 5 . The reasons for 
this enhancement also apply here. 

It is well known that the strength of strong penguins is an important issue of current interest, 
relevant to CP violation, tests of the SM and to searches for new physics signals at BABAR and 
BELLE. The first evidence for a strong penguin was found by CLEO in 1997 by measuring B — > 



Krj and B — > r]'X s [39], which turned out be "unexpectedly large". Recently both BABAR H] 
and BELLE [FJ have confirmed the CLEO measurement with improved precision. Theoretically, 
there are many models to explain the yield of rj' in B decays. The measurement of the rj' 
momentum spectrum by CLEO |3(J has ruled out models with a b — > ccs enhancement through 
a possible cc content in the r( wave function. Up to now, there are two surviving explanations 
for the high yield of rj in B decays: A model in which new physics enhances the strong penguin 
b — ► sg and a model incorporating the QCD anomaly coupling g — g — r! [^, |33|]. The rj is 



a very complicated object in QCD. There are unsolved puzzles involving the ?/, such as its large 
mass, mixing between flavor singlet and octet, gluonic content, decay constants and its mixing 
with glueballs, which inhibit reliable theoretical predictions for B — > rjK^*' and B — > vjX s . 
The aforementioned difficulties are absent for <fi and theoretical predictions for B — > <pX s are 
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rather clean. Furthermore, <\> has a clear experimental signature. We can therefore conclude 
that the theoretical study and experimental measurement of B — > <pX s will provide quite a 
clean ground for testing the SM and searching for new physics. 

In summary, we have studied the semi-inclusive decays B — > <pX s . We have calculated 
the contributions of b — > <fts and of the b — > <psg subprocesses with the gluon radiated from 
the charm loop or from the off-shell gluon of b — > sg*. These effects are found to be quite 
large, giving B{b — > <psg) = 3.8 x 10~ 5 . Cutting the <p momentum |k^| > 2.0 GeV, results in 
B{b — > (psg) = 1.0 x 10~ 5 . Adding this contribution to the dominant process b — > (f)s, we predict 

B(B -> <j)X s ) = 10.5 x 10" 5 (29) 

and 

B(B -> 4>X S ) = 7.1 x 10" 5 (|k | > 2.0 GeV). (30) 

This large decay rate and its clear signature render detailed studies of B — > <j)X s at BABAR and 
BELLE feasible in the near future. These will shed light on the strength of strong penguins. 
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Appendix A: The operators Oi 

The operators in Eq.[l| read 





= c a ^Lb a ■ sp^^Lcp , 


o 2 


= Cofl^Lbp ■ 




o 3 


= s a ^Lb a ■ spj^Lsp , 


o A 


= Sa^Lbp ■ 




o 5 


= Sa^Lba ■ , 


o 6 


= SaYLbp ■ 




o 7 


= \s a ^Lb a ■ esS^^Rsp , 


o 8 


= paYLbp 


■ e s s^^Rs a 


o 9 


= \s a ^Lb a ■ egSpry^Lsp , 


Oio 


= paYLbp 




o g 


= ^s a a^RT^m b b p G%. 









(31) 
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Where a and (3 are the S77(3) color indices and L — (1 — 7s)/2, i? = (1 + 75)/2. 



Appendix B: Functions for the charm loop 

Using the notations ri = r 2 — q for the momentum of the on-shell gluon in the final 

state and k for the momentum of the off-shell gluon, we have 

Az 5 = -2 + 16 + 2r2 (G (r 2 ) - G (n + r 2 )) - - (G_(r 2 ) - G-(ri + r 2 )) , (32) 
ri ri 



7 7 r 2 r\ r\ 4r 2 (l + r 2 - r^) 

A ^6 = J + — +o 2 G 0(^2) 

4 ri 2ri 3rf 3ri rf 

2(r 1 + 3r* + r 2 + r 1 r 2 ) r , 4(r x -3r 2 ) ^ 4(1-372) , , , 

H 2 G (ri + r 2 ) 5 G -( r 2) H G -( r i + r 2) 

rf rf r\ 

+ 5(4-^ r + 5(r? + 2 ri (r, - 2) + r 2 (r 2 - 4)) + 

rf rf 

4 2r 
A* 23 = -2 + -(G_(r 2 )-G_(r 1 +r 2 ))-^(G (r 2 )-G (r 1 + r 2 )), (34) 

A* 25 = -2(G (r 2 )-G (r 1 +r 2 )), (35) 

8, H 22 16 + 2ro . . 16-4n + 2r 2 _. . 

Az 2 = - ln-^ + — + — G (r 2 ) G (n + r 2 ) 

3 m c 9 3ri 3ri 

--(G_(r 2 )-G_(n + r 2 )), (36) 
n 

A* 85 19-3r 2 r2 /l \ 4r 2 / 3 \ _ , N 

A* 3 = 41n-^-- + * + -*- --1 +-1 l-r 2 -- ) G (r 2 ) 

m c 3d ri 9ri \r\ J Sri \ t\) 

/10 2 14r 9 2r 9 4r?\ „ , 4 / 2r 9 



^2 ^2 _ **' 2 
2 + ^ _ Q^2 



l) G (r 1 + r 2 ) + -(l + G_(r 2 ) 



\ 3 3ri 3rf r\ 3r 

-- (l + — - r 2 ) G?_(n + r 2 ) + ^(4 - r 2 )T (r 2 ) 
ri \ ri / rf 

+9 ( 1 _±_^ + ^ + lfW (ri + r2 ), (37) 



Ai 8 = H ]n iL_^ + 8 (r2 _4)G„(r a )-|(l-- + ^)G„(r I+ r 2 ), (38) 
3 m c 9 3ri 3 V ri ri/ 

A in = -| In -g- + H + l Go ( r2 ) + 2(8 + J' + r2) (Gb W - G (n + r 2 )) 
6 m c 9 6 3ri 

(2G_(n) - G_(r 2 ) - GL(n + r 2 )) , (39) 

Az' 12 = In A + 4 + 24 - 2r2 - ^ [(1 + 2r!)r 2 G (r 2 ) - (rf + r 2 + 2r 1 r 2 )G (r 1 + r 2 ) 

3 m c rx rf L 

+ ^(G_(r 2 )-G_(r 1 + r 2 )) 
rf 
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+- 



12 r 



r 2 (r 2 - 4)T (r 2 ) + (rf + 2n(r 2 - 2) + r 2 (r 2 - 4))T (n + r 2 ) , (40) 



A* 15 = -\ (l + -)G (r 2 ) + Ul + ^— )G ( ri + r 2 ), (41) 
3 V r 2 / 3 V ri + r 2 / 

2 2 (8 + r 2 )^ 2 .8 + r 2 4 

Az i7 = -Go(r 2 ) + -( + — — )G (ri + r 2 ) 

3 3ri 3 ri ri + r 2 

+^(G_(r 2 )-G_(r 1 + r 2 )), (42) 

A* 19 = -^ + ^(3 + ^--)(G (r 2 )-G (ri+r 2 ))+ 1 , G (ri + r 2 ) 
3 3 ri ri 3(ri + r 2 J 

+ ^(G_(r 2 )-G_(r 1 + r 2 )), (43) 



A . _ 2 16 8 r 2 2r 2 (24n + 3r| + 20r 2 + lr x r 2 + 4r^ g 
3 r\ 3 ri 3r^(ri + r 2 ) 

+ ^|(20 + 3n + 4r 2 )G (r 2 ) + - 2 [(n + 4r 2 )(G_(r 1 + r 2 ) - G_(r 2 )) 
+2r 2 (4 - r 2 )T (r 2 ) + 2(n + r 2 )(n + r 2 - 4)T (n + r 2 )] . (44) 



The function Gi(t)(i = —1, 0, 1, 2) are defined as 

GAt) = f 1 dxx i lnfl - tx(l - x) - iS\. (45) 
Jo 

Gi )2 (t) can be simplified to Go. The explicit form of G-i t o(t) could be found in Ref. Jl9|. In 
the equations above we used G-{t) = G-\(t). 
To,i are defined by 



x 



Ti = / dx - r (46) 

Jo 1 — x(l — x)t — ie 



Their explicit forms are given by 



4arctan V — ■ 0<t<4 



T (t) = { r ' (47) 

2 i7r+21n(v / 7-v / ^4)- 2 lri(v / T+y / F : 4) . t > A 

y/t{t-4) 

The one loop two-point charm penguin function is well known in the form of 

G(t) = -\ + ~ In JL - f 1 d^^l - ln[l " £(1 " Ot + M], (48) 
3 3 m c Jo 

whose explicit form is given by 

G(t) = - 8 - + t^^-l( l + l)Co(t). (49) 



9 3 m r 3 V i 
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Appendix C: Squares and interferences of amplitudes 



Now we write down the terms in Eq.|25] needed for calculating the inclusive decays 



\M f 



\M f 



57T 



f \Cs\ 2 Sml 



4| J F 1 | 2 st 2 (s + t) + 2$(F F*)tu(-7s 2 + 6st + 5t 2 



+ \F \ 2 u(l2st(t + u)+ 4s 2 (7t + u) + t 2 {t + llu)) 



(50) 



|F 3 | 2 (s 3 + 4t 2 w + s 2 (w - 2t) + st(* + 5u 



-2^(F 3 F;)s(s-t) + \F 5 \ 2 s(s + t) 



(51) 



I7W 



For the charm loop penguin contributions, we get 

2 / r \ 2 

|C 1 | 2 2m 2 [4s 2 (s + u)\Jl\ 2 + 4t(tu + s(t + 2u))\J\ 



49 / <? 



A 



81 V 16tt 2 



4 



|2 
i6l 



+tw 2 (2st + u(s + t))\ J° 23 | 2 + t(2s 3 - st 2 - t 3 )| J/ 23 | 2 



+5? 



rl* 
'i23 



4J°(2s 2 (s + t)j£ + 2s(tu + at + su)J°* + stu{s + u) 4* + s 3 wJ/ 2 * 
- S 2 t( S + t)4 3 ) + 4*4(2t(a + t)4* - (s 2 - - *>4* 3 + *t( a + t) J* 
+4J°( M (t 2 M - s{t + u)(s- 0)4*3 + 2s 2 tuJ% 3 - st 2 (s + i) 4* 3 ) 
+2suJ° 23 (s(tu -st- su)Jl* z + t{s 2 -st- t 2 ) J 2 * 3 ) - 2s 2 t 2 (s + *)4 3 4* 3 ] } (52) 

(t^) 2 (t) 2 '^' 24 ^ 2 



X 



2s 2 ( s + «)|4| 2 + 2st(t + «)|4| 2 + 2tu{s + t)i4 2 | 2 + Asm 044* - uJgj\ 



'il2 



-4s 2 (s + t)^44* - 2s 2 twK 



tO t1* 
J i8 J il7 



TO t1* I Tl TU* Tl T' 



1 tO* 



rl tO* 



-stu^ft 



2(s + «)4 2 4 7 + 2t4 2 4*i + *(* + «)4r4r 



tO* 



'i%°GX i8°il7 



rl tO* 



i8^«21 



s t JnjJi2\ + JinJi2i 



The interference terms are read as 



(53) 



K(M a+ .mL 



|dC 8 | 2 ^ 8m 2 ^|2 S t[F 3 ( S -t)- F 5 (s + 



-F*[8F 4 s 2 (s + t) + F 3 u(2s 2 - 9st - 3t 2 ) + F 5 u(3st + t 2 - 2s 2 )]|, (54) 



(5 



2s 2 (4 + 4) + 2^(4* + 4) + 2n( s + 1)4 + stu{Ji* 7 - j t 



rO* n 



rO* 



0* 



rl* ~> 

«21y 
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$t(M A+ MU 



U(M A+ Ml 



to(M A _MU) = 



2(5 + t)J% + tu(J% 7 - + 44 [si 2 4 + (st(s + t) - s 2 u)4 



st(s + u) 4] + 2st«4 + (t + «) Jgi + t{J& - 4* 
+mJ° 6 [2s(s 2 - st - t 2 ) J}* - 2t(s 2 - su + fu) j£ 2 - st ( s2 + tM M 
-2st 2 4* - st(st + su- ut)4* - 2st{s + u)J\ 



i* 

(17 



tO* 
>i2 



n 

J i26 



2s 2 t( s + t)Jl; 



/o* 
£17 



-2su(s 2 + st + t 2 ) J°* + 2t 2 w(2s + t) 4* 2 + stw(s 2 - iu) 
+s 2 t 2 u(J^ 7 - j} 2 \) + stu(tu + st + sit) j, 

/ _ v 9 

S 2 + 2tM + SM - St)(sJ° 5 + t J° 6 ) 



(55) 



4m?K —F, 



9 128vr 4 1 V 4 7 

+2st M (sJ° 3 + t4 3 ) + 4 6 (s(t + m) - s 2 - 2*u) - 2t\(s + t)4 
\s + t)(2sJ° + 2tJ° + stjig) - u(s 2 J° 3 + (st + t 2 ) J° 6 )^ 
7° " 

y i26 



;26 



-2F 4 *s 



+F 5 *s(2s J° + 2*4 - tu 4, 



(56) 



7 fl2 



9 128tt 



: |CiC 8 | ^ 4m 2 ^{F;st(4sJ° + 4* 4 + 3snJ° 3 - 2tnJ° 6 



-F * [4s 2 (s + t)(tJ 2 23 - 24) - s 2 u(4s + t)4 3 + 4us(s - t)J t 
+ (4s 2 + 6st - t 2 )tuJf 2e - 2(4s 2 - st + t 2 )uJi 6 + stw 2 4 
+uH 2 (t-4s)J° 2 



o 

'a 



(57) 



1287T 



t M 4i 



s 2 (s + t)4 - 2tn(5s + 2t)4i2 - stu(3s + t)4i 7 + tu(As 2 + st - t 2 )J} 21 

] (58) 

s 2 (t-s-n)44 
s + M )4 7 ) 



128tT 



F x *st[4s4 + 4*4 - 3suJ° 7 - 3^4^ } 
0} Zmfrfa 



tu(2(s + t)J? 12 -s(2J° 8 + tJ? 21 + 

2(s + t)4 + M4 7 - j; 



st(t + u- s)4 3 



-F 4 V 



■o - 

i21, 



F*sfs4 + t4 



(59) 



The functions Fi are defined by 



F 3 



x)ml 



o ^ kfik 2 ' 
1 x(p(x)mf 







A' 2 



(60) 
(61) 
(62) 
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Fi 

Fr, 



o 



dx- 



/c 2 
(x)mj 



G(r 8 ), 



(63) 
(64) 



and r3 



fx) is the distribution function of 6 meson, which can be found in Ref. 28 . In 



numerical calculation we take 4>(x) = 6x(l — x), which kills the IR pole of the gluon propagator 
p. The other gluon is always time-like k$ > 0, so there is no IR divergence in our calculations. 
The J functions are 



7° 

J i2 



7° 



7° 



dxTn 



dxTn 



(x)Ai 2 



r 2 
{x)Ai 5 



dxro 



r 2 
[x)Ai 6 



7° 

J il2 



dxr. 



r 2 

xcf)(x)Ai s 



dxrc 



r 2 
{x)Ai 12 



o 



r 2 



J il7 



J i21 



n 



o 



1 , 2 2x(fi(x)Ai 17 
dxr 



dxr$ 



dxr$ 



r\r 2 
2x(p{x)Ai 21 



r\r 2 
2x(f)(x)Ai 23 



7° 

J i26 



dxr\ 



r x r 2 
2<f){x)Ai 26 



r\r 2 



7° 



J *=/„ 

Jo 

T° - 



7° 

J ill 



7° 

J il7 



J 



i21 



7° 

^23 



7 2 

^23 



/2(i 



dxr n 



dxri 



dxvQ- 



dxrn 



[x)At 3 

?~2 

x0(x)A2 5 
(x)Az 8 
J(z)Ai u 



dxr 2 . 



^2 

2</>(x)Aii 7 



dxr, 



rir 2 

2 2(j)(x)Ai 21 

r\r 2 

2 20(x)Ai 23 
o 



rir 2 

i 2 2x 2 0(x)Ai 23 
ctxrn 











r x r 2 



1 2 2x0(x)Ai 26 
rir 2 



(65) 
(66) 
(67) 
(68) 
(69) 
(70) 
(71) 
(72) 
(73) 
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